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ABSTRACT 


Measured,  local-average  Nusselt  numbers  (Nu)  for  in-tube,  convective  boiling  and  condensation  with  a 
twisted-tape  insert  are  presented  for:  R12,  R22,  R152a,  R134a,  R290,  R32/R134a,  R32/R152a, 
R290/R134a,  and  R134a/R600a.  The  heat  transfer  data  were  obtained  from  a fluid  heated/fluid  cooled, 
9.7  m long  condenser  and  evaporator  of  a breadboard  refrigeration  cycle.  Convective-boiling,  heat- 
transfer  data  were  taken  for  transition  and  turbulent  all-liquid  Reynolds  numbers.  Convective- 
condensation,  heat-transfer  data  were  taken  for  laminar  and  turbulent  all-liquid  Reynolds  numbers.  The 
measured  convective  boiling  and  condensation  Nusselt  numbers  for  the  single  component  and  the 
azeotropic  mixtures  were  each  correlated  to  a single  expression  consisting  of  a product  of  dimensionless 
properties.  The  single  component  convective-boiling  correlation  was  modified  to  predict  the  zeotropic 
mixture  data.  The  predictions  obtained  from  the  modified  flow-boiling  correlations  found  in  the  literature 
were  significantly  different  from  the  present  Nusselt  number  measurements.  Presumably,  the  correlations 
from  the  literature  could  account  for  neither  the  partial  dryout  induce  by  the  tape  nor  the  bubbly-mist  flow 
introduced  by  an  expansion  valve. 


Keywords : Enhanced  heat  transfer;  refrigerant  mixtures;  flow  boiling;  in-tube  condensation;  twisted  tapes; 
building  and  technology 
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NOMENCLATURE 


English  symbols 

a constant  in  eqns.  1 and  2 (K) 

Ac  cross-sectional  flow  area  (m2) 

b constant  in  eqn.  2 (K/m) 

a" 

Bo  local  boiling  number,  — A — 

Gnfg 

c constant  in  eqns.  1 and  2 (K/m2) 

cp  specific  heat  (J/kg»K) 

Dh  hydraulic  diameter  (m) 

Dj  inner  diameter  of  tube  (m) 

Eh  relative  U (%)  in  h 

f Fanning  friction  factor 

G mass  velocity  of  empty  tube  (kg/m2  • s) 
h local-averaged  heat-transfer  coefficient  (W/m2  • K) 
hfg  latent  heat  of  vaporization  (kJ/kg) 

H 180°  twist  pitch  (m) 

Ja  Jakob  number,  cp,»ATs/hfg 

k liquid  thermal  conductivity  (W/m  • K) 

Nu  local-average  Nusselt  number  based  on  D, 
m mass  flow  rate  (kg/s) 

M gram-molecular  weight  (g/mole) 

P local  fluid  pressure  (Pa) 

Pc  critical  pressure  (Pa) 

Pr  local  liquid  Prandtl  number 

q"  local  heat  flux  (W/m2) 

P 1 GD  : 

ReD  all  liquid,  empty  tube  Reynolds  number,  — — - 

Pi 

Res  all  liquid,  swirl  Reynolds  number 

Sw  swirl  parameter 

T temperature  (K) 

UA  overall  conductance  (W/K) 

U expanded  uncertainty 

itl 

Va  mean  axial  velocity,  (m/s) 

Pr^c 

Vs  swirl  velocity  (m/s) 

x mole  fraction 

xq  thermodynamic  mass  quality 

X normalized  axial  heat  exchanger  coordinate 

y twist  ratio,  H/D; 


Greek  symbols 

5t  thickness  of  twisted  tape  (m) 

AL  heat  transfer  length  of  one  increment  (m) 

ATS  wall  superheat  or  wall  subcooling,  | Tw  - Ts  | (K) 

Aq  duty  of  one  element  (W) 

0 dimensionless  mixture  temperature  difference,  (Td  - Tb)/(TLV  - Tmv) 
/x  liquid  dynamic  viscosity  (kg/m  • s) 
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X, 


Martinelli  parameter 


subscripts 

a annulus 

b bubble  point 

c condensation 

d dew  point 

e boiling 

1 liquid,  liquid  composition 

LV  less  volatile  component 

m mixture 

MV  more  volatile  component 

p single  component\azeotrope 

r parameter  normalized  by  critical  value 

s saturated  state 

t twisted  tape 

v vapor,  vapor  composition 

w inner  tube  wall 

x local-averaged  quantity 
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INTRODUCTION 

Twisted-tape  inserts  have  been  used  to  enhance  heat  transfer  since  the  19th  century.  Marine  steam  boilers 
were  fitted  with  "retarders"  (twisted  tapes)  to  reduce  coal  consumption.  Research  began  to  appear  in  the 
literature  around  this  time.  In  one  of  the  first  studies,  Whitham  (1896)  showed  that  twisted-tapes  resulted 
in  a 18.4%  fuel  savings  for  a coal-fired,  horizontal,  tube  boiler  of  a Philadelphia  railway  station.  In  a 
more  recent  study,  Mamer  and  Bergles  (1978)  reported  300%  enhancements  over  the  heat  transfer  of 
empty  tubes  for  single  phase  flow  with  twisted-tape  inserts.  In  the  time  between  the  Whitham  (1896)  and 
the  Mamer  and  Bergles  (1978)  publications,  90  manuscripts  were  published  on  the  enhancement  of  single 
phase  convection  with  swirl  flow  devices  (Bergles  et  al.,  1995).  As  of  1995,  the  number  of  publications 
on  the  subject  had  grown  to  261.  Over  the  years,  the  vibrant  interest  in  the  study  of  twisted  tapes  as  a 
heat  transfer  augmentation  has  been  sustained  by  its  affordability  and  its  suitability  for  retrofiting  existing 
shell-and-tube  heat  exchangers. 

The  flow  enhancement  of  the  twisted  tape  arises  primarily  from  increased  flow  path  length  and  swirl 
mixing.  Swirl  flow  and  increased  path  length  are  expected  to  also  benefit  two-phase  heat  transfer.  For 
example.  Royal  and  Bergles  (1976)  surveyed  horizontal,  convective  condensation  with  twisted  tapes  and 
found  improvements  in  heat-transfer  coefficients  by  as  much  as  30%  over  empty  tube  condensation.  In 
addition,  the  most  popular  use  of  twisted  tapes  in  flow  boiling  is  to  delay  the  occurrence  of  burnout 
(Gambill  et  al.  1961,  Bergles  et  al.  1971,  and  Papadopoulos  et  al.  1994.)  Yet,  relatively  few  two-phase 
heat  transfer  studies  with  twisted  tapes  exist:  approximately  ten  convective  condensation  and  possibly  60 
flow  boiling.  Furthermore,  to  the  best  of  the  authors’  knowledge,  no  two-phase  heat  transfer  studies  with 
alternative  refrigerants  nor  refrigerant  mixtures  for  twisted  tapes  exist.  This  is  unfortunate  considering 
that  using  twisted  tapes  is  an  inexpensive  way  to  reduce  heat  exchanger  size  or  improve  cycle  efficiencies. 
Consequently,  there  is  a need  for  phase  change  heat  transfer  data  that  can  be  used  to  design  heat 
exchangers  with  twisted-tape  inserts  for  alternative  refrigerants  and  mixtures. 

This  paper  presents  measured  local-averaged  heat-transfer  coefficients  for  in-tube  convective  boiling  and 
condensation  of  several  possible  alternative  refrigerants  with  a twisted-tape  insert.  The  heat  transfer 
performance  of  five  single  component  and  two  azeotropic  refrigerants  were  examined:  R12,  R22,  R152a, 
R134a,  R290,  R290/R134a  (0.44/0.56  % mass),  and  R134a/R600a  (0.81/0.19  % mass).  Also,  the 
performance  of  two  zeotropic  refrigerant  mixtures  - R32/R134a,  and  R32/R152a  - were  examined  at 
approximately  five  different  mass  compositioas  of  R32.  Table  1 presents  the  various  percent  mass 
compositions  of  R32  at  which  the  zeotropes  were  tested.  The  heat  transfer  tests  of  the  various  working 
fluids  were  used  in  an  attempt  to  develop  universal-fluid  evaporative  and  condensing  heat  transfer 
correlations. 

Improved  evaporator  and  condenser  designs  can  be  realized  with  the  use  of  local  heat  transfer 
correlations.  Local  heat  transfer  correlations  can  be  used  to  tailor  evaporators  and  condensers  for  quality 
ranges  and  flow  rates  that  are  germane  to  a particular  application.  True  local  heat  transfer  measurements 
could  have  been  obtained  with  an  electrically  heated  test  section.  However,  electric  resistance  heating 
is  not  a physically  realistic  boundary  condition  for  refrigerant  cycle  applications.  Accordingly,  local- 
averaged  heat-transfer  coefficient  measurements  were  obtained  from  a fluid  heated  test  section  that  was 
divided  into  several  segments.  The  heat-transfer  coefficient  was  averaged  over  a particular  range  of 
qualities  for  each  segment  providing  quasi-local  (local-averaged)  data.  Kattan,  et  al.  (1995),  Goto,  et  al. 
(1995),  and  Conklin,  et  al.  (1992)  provide  examples  of  local-averaged  flow-boiling  experiments.  An 
additional  uncertainty  in  the  heat-transfer  coefficient  is  introduced  by  the  averaging  process.  If  the 
variation  of  the  heat-transfer  coefficient  over  the  quality  range  for  which  it  was  averaged  is  sufficiently 
small,  the  uncertainty  is  negligible. 

TEST  APPARATUS 

Figure  1 is  a schematic  of  the  breadboard  refrigeration  cycle  from  which  the  heat  transfer  measurements 
were  made.  The  figure  shows  the  location  of  the  evaporator  and  condenser.  The  evaporator  and 
condenser  each  consisted  of  twenty  0.5  m lengths  of  copper  spined  pipes  with  a snugly  fitting  copper  tube 
over  the  spined  fins.  The  twenty  horizontal  segments  were  stacked  vertically  and  joined  in  series  by  u- 
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bends  and  T-connections  as  shown  in  Fig.  2.  The  test  refrigerant  flowed  through  the  9.64  mm  inner 
diameter  tube  which  contained  an  aluminum  twisted-tape  insert.  The  180°  twist-pitch-to-tube-diameter 
ratio  (twist  ratio,  y)  was  4. 15  and  the  tape  thickness  (6)  was  0.5  mm.  A 40%  mass  glycol/water  mixture 
flowed  through  the  spined  annulus  counterflow  to  the  refrigerant.  The  refrigerant  was  delivered  to  the 
entrance  of  the  evaporator  by  a thermostatic  expansion  valve  at  approximately  12%  quality.  Refrigerant 
vapor  exited  the  evaporator  with  approximately  14  K of  superheat.  An  open-drive,  reciprocating 
compressor  supplied  superheated  vapor  to  the  condenser  entrance.  The  refrigerant  exited  the  condenser 
as  subcooled  liquid.  Although  neither  an  oil  separator  nor  an  accumulator  were  used  in  the  cycle,  the 
measured  oil  concentration  was  below  the  accuracy  of  our  measurement  which  was  0.1%.  Further  details 
of  the  cycle  can  be  found  in  Pannock  and  Didion  (1991.) 

Figure  2 illustrates  the  flow  paths  in  the  heat  exchanger  and  the  location  of  the  refrigerant  and  the 
secondary  fluid  temperature  measurements.  The  geometry  and  the  size  of  the  evaporator  and  condenser 
were  practically  identical.  Thermocouples  were  epoxied  to  the  outside  of  the  refrigerant  tube  bends  and 
insulated.  Instream  thermocouples  were  located  in  ever  other  bend  of  the  secondary  fluid  stream.  A ten- 
element,  instream  thermopile  was  used  to  measure  the  temperature  change  in  the  secondary  heat  transfer 
fluid  across  the  entire  heat  exchanger. 

Fluid  mixers  were  not  used  before  the  thermopile  ends.  From  past  experience  when  mixers  were  used 
in  other  experiments,  the  expanded  uncertainty  (U)  of  the  temperature  difference  measurement  with 
thermopiles  was  estimated  to  be  within  0.001  K to  0.01  K.  (The  expanded  uncertainty  estimates  were 
obtained  from  a root-sum-square  of  the  component  uncertainties  with  a 95%  confidence  interval.)  Two 
bends  and  0.2  m of  unheated/insulated  length  at  the  exit  of  the  heat  exchanger  provided  some  mixing  and 
time  to  approach  a uniform  temperature  before  the  thermopile.  The  inlet  fluid  was  mixed  well. 
Considering  that  the  mixed  state  of  the  heat  transfer  fluid  at  the  heat  exchanger  exit  was  not  as  certain 
as  the  inlet,  the  expanded  uncertainty  (U)  was  estimated  to  be  0. 1 K. 

The  conduction  and  mixing  errors  introduced  by  the  indirect  measurement  of  the  refrigerant  temperature 
at  the  insulated  tube  bend  wall  can  be  relatively  small  for  a thin  tube  wall  and  small  driving  temperature 
differences.  Errors  in  the  refrigerant  temperature  measurement  due  to  conduction  of  heat  along  the  tube 
wall  from  the  active  portion  of  the  heat  exchanger  to  the  insulated  tube  bend  were  modeled  as  one- 
dimensional heat  conduction.  The  conduction  error  estimated  using  this  method  typically  ranged 
(depending  on  the  q")  from  0.0009ATS  (=  | Tw  - Ts  | ) to  0.009 ATS  or  within  0.3  K for  most  heat  fluxes. 
Mixing  errors,  those  caused  by  a nonequilibrium  refrigerant  state  at  the  point  of  temperature 
measurement,  are  due  to  insufficient  mixing  of  the  superheated  sublayer  in  the  adiabatic  tube  bend.  Swirl 
and  bend  flow  should  contribute  to  some  mixing  of  the  refrigerant  to  an  equilibrium  state.  However,  a 
maximum  mixing  error  can  be  estimated  if  it  is  assumed  that  no  mixing  occurs  in  the  bend,  and  the 
problem  is  simplified  to  a transient  conduction  problem.  Using  this  method  the  liquid  sublayer  was 
estimated  to  come  within  0.08ATS  of  the  saturation  temperature  by  the  time  it  reached  the  apex  of  the  tube 
bend.  Considering  that  ATS  varied  on  average  from  1 K to  11  K for  boiling  and  from  1 K to  6.5  K for 
condensation,  the  maximum  error  was  estimated  to  be  within  0.88  K and  0.52  K for  boiling  and 
condensation,  respectively.  Combining  the  conduction  and  mixing  errors,  the  maximum  error  was 
estimated  to  be  within  0.92  K.  Considering  that  the  swirl  and  bend  flows  contribute  some  mixing,  the 
maximum  error  should  be  less  than  0.92  K. 

The  mass  flow  rates  of  the  40%  mass  glycol/water  mixtures  were  measured  with  Coriolis  flow  meters 
having  an  expanded  uncertainty  of  0.001  kg/s.  Turbine  flow  meters  were  used  as  a secondary  flow  rate 
measurement.  The  turbine  meter  measurements  agreed  with  the  measurements  obtained  from  the  Coriolis 
meters  to  within  the  uncertainty  of  the  turbine  meter.  The  refrigerant  mass  flow  rate  (m,)  was  calculated 
from  an  energy  balance  on  the  entire  evaporator.  The  relative  expanded  uncertainty  of  rii,  was  2 % to 
3 % . The  average  thermodynamic  quality  (x,,)  for  each  heat  exchanger  element  was  calculated  from  an 
energy  balance  on  that  element.  The  relative  expanded  uncertainty  of  Xq  was  within  7%  of  the 
measurement.  Table  2 shows  the  estimated  relative  uncertainty,  U (%),  for  various  measurements. 
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DATA  REDUCTION  METHODOLOGY 

Figure  3 shows  the  refrigerant  and  secondary  fluid  temperature  profiles  for  a R22  evaporator  and 
condenser.  The  temperatures  are  plotted  against  a dimensionless  coordinate  (X)  which  represents  a 
fraction  of  the  axial  length  of  the  heat  exchanger.  For  example,  the  refrigerant  entrance  to  the  evaporator 
and  the  condenser  is  shown  at  X = 0 and  X = 1,  respectively.  At  approximately  X = 0.88,  saturated 
refrigerant  vapor  exists  for  both  the  evaporator  and  the  condenser.  The  vapor  is  superheated  for  values 
of  X greater  than  approximately  0.88.  Heat-transfer  data  reported  here  are  only  for  the  two-phase  region, 
i.e. , for  thermodynamic  qualities  between  0 and  1 . 

The  uncertainty  in  the  Nusselt  number  was  reduced  by  using  a regression  of  the  secondary  fluid 
temperatures  to  X in  its  calculation.  Figure  3 shows  that  the  scatter  in  the  evaporator  secondary  fluid 
temperature  profile  is  marginally  greater  than  that  exhibited  in  the  condenser  secondary  fluid  temperature. 
A weighted  regression  was  used  to  reduce  the  uncertainty  caused  by  random  effects  in  secondary 
temperature  measurement.  Those  evaporator  secondary  fluid  temperatures  not  associated  with  the 
superheated  refrigerant  vapor  were  fitted  to: 

Tf  = a + cX2  (1) 

The  weighted  regression  takes  advantage  of  the  stability  of  the  inlet  secondary  fluid  temperature 
measurement.  Unfortunately,  the  inlet  secondary  fluid  temperature  was  omitted  from  the  fit  because  it 
lies  in  the  superheated  refrigerant  region.  To  overcome  this,  the  exit  secondary  fluid  temperature  was 
calculated  from  the  inlet  temperature  and  the  temperature  change  obtained  from  the  thermopile.  The 
variance  of  the  calculated  exit  fluid  temperature  was  approximately  25  times  smaller  than  that  of  the 
secondary  fluid  temperature  measurements  taken  in  the  heat  exchanger  bends.  Accordingly,  a weighing 
factor  of  25  was  used  on  the  exit  fluid  temperature.  Only  a 0.3%  change  in  the  model  coefficients 
occurred  when  the  weighing  factor  was  decreased  from  25  to  15.  Consequently,  a negligible  error  is 
introduced  for  a large  error  in  the  estimate  of  the  weighing  factor. 

A weighing  factor  was  not  used  for  the  condenser  calculation  because  of  the  relatively  small  scatter  of 
the  secondary  fluid  temperature  measurements  and  because  the  exit  and  inlet  secondary  fluid  temperatures 
were  not  associated  with  two-phase  refrigerant  flow.  Nevertheless,  the  secondary  fluid  temperatures  were 
fitted  to  X to  reduce  the  uncertainty  in  the  measurement.  As  shown  in  Fig.  3,  the  form  of  the  condenser 
secondary  fluid  temperature  differs  from  that  of  the  evaporator.  Consequently,  the  secondary  fluid 
temperatures  associated  with  the  two-phase  refrigerant  of  the  condenser  were  fitted  to: 

Tf  = a + bX  + cX2  (2) 


The  forms  of  eqns.  1 and  2 provided  the  best  fit  of  the  measurements  with  the  minimum  number  of 
terms. 

The  overall  conductance  (UA)  equation  can  be  written  for  a single  0.5  m long  heat  exchanger  element 

as: 


UA  Aq  htAt  2tc  ALkcu  T)haA0 


The  h,  is  the  heat-transfer  coefficient  of  the  inner  smooth  tube  with  the  twisted-tape  insert  and  is  based 
on  the  empty  tube-side  area  A,  (xDjAL.)  The  next  to  the  last  term  is  the  thermal  resistance  of  the  tube 
wall  between  the  fluids.  The  h,  represents  the  heat-transfer  coefficient  for  the  spine-fin  annulus. 
Kedzierski  and  Kim  (1996)  developed  a single  phase  heat  transfer  correlation  for  a spine-fin  annulus  with 
the  same  geometry  as  the  spine-fin  annulus  in  the  present  test  heat  exchangers.  It  was  presumed  that  the 
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bends  destroyed  the  fully  developed  thermal  boundary  layer.  Consequently,  the  1^  correlation  was  used 
to  account  for  thermally  developing  flow  in  each  element. 


The  log-mean  temperature  difference  (AT^  was  calculated  for  each  increment  using  the  measured 
refrigerant  temperatures  and  the  secondary  fluid  temperatures  from  eqn.  1 . The  effective  heat  capacity 
of  a zeotropic  mixture  is  nearly  constant  for  each  increment.  Consequently,  the  error  introduced  in  using 
the  log-mean  temperature  difference  equation  was  negligible. 

The  duty  for  the  ith  element  was  calculated  from  eqn.  2 as: 


(4) 


Aq  = ifi£  CPf  (bCXj.,  - Xj)  + c(Xi<£2  - X£2)) 


where  X,  and  Xi+1  are  evaluated  at  the  refrigerant  inlet  and  exit  of  each  increment  of  length  AL,  and  b 
is  zero  for  the  evaporator.  All  of  the  core  elements  were  0.5  m in  length.  The  length  of  the  end 
elements  were  shortened  in  an  appropriate  amount  to  account  for  the  single  phase  heat  transfer  in  the 
element. 

The  local-average  heat-transfer  coefficient  (h,)  was  calculated  for  eqn.  3 as: 


/ 


(5) 


The  average  expanded  uncertainty  for  the  h,  was  within  9%  for  all  of  the  data. 

TEST  RESULTS 

Appendix  A contains  four  tables  of  all  the  measured  data  generated  in  this  study.  The  azeotropic/pure 
flow  boiling  and  convective  condensation  data  are  given  in  Tables  A1  and  A2,  respectively.  The 
zeotropic  flow  boiling  and  convective  condensation  data  are  given  in  Tables  A3  and  A4,  respectively. 
All  the  parameters  are  defined  in  the  nomenclature.  The  friction  factor  (f),  and  the  Martinelli  parameter 
(Xtt.)  were  not  used  in  the  correlation  of  the  data,  but  are  given  in  the  spirit  of  completeness.  The  Ej,  is 
the  relative  U in  the  heat-transfer  coefficient.  Data  with  values  of  Eh  greater  than  0.125  were  omitted 
from  the  correlation. 

Over  200  graphs  would  be  required  to  plot  the  measured  local-averaged  heat-transfer  coefficient  (h,) 
versus  thermodynamic  quality  for  each  test  run.  Instead,  only  representatives  plots  of  1\  are  presented 
in  Figs.  4 through  21  for  each  fluid  under  flow  boiling  and  convective  condensation  conditions.  Several 
different  predictive  correlations  are  compared  to  the  measured  data.  A more  detailed  discussion  of 
predictions  and  measurements  follows. 

Convective  Condensation 

Figures  4 through  12  show  the  measured  local-averaged,  convective-condensation,  heat-transfer  coefficient 
(h,)  as  a function  of  thermodynamic  quality.  The  data  are  compared  to  predictions  from  four  different 
correlations.  No  universal-fluid  local  convective-condensation  coefficients  for  twisted  tapes  were  found 
in  the  literature.  Three  universal-fluid,  smooth-tube,  horizontal,  convective-condensation  correlations 
were  modified  according  to  a procedure  given  by  Royal  and  Bergles  (1978.)  The  procedure,  as  explained 
below,  permits  the  smooth-tube  condensation  correlations  of  Soliman  (1968),  Shah  (1979),  and  Traviss 
et  al.  (1973)  to  account  for  the  enhancement  of  the  twisted  tape.  The  fourth  correlation,  shown  as  a solid 
line,  is  a fit  of  the  measurements.  No  general  convective  condensation  correlations  for 
azeotropic/zeotropic  mixtures  were  found  in  the  literature.  Consequently,  the  mixture  properties  were 
used  in  the  pure  component  condensation  correlations  to  predict  the  mixture  heat  transfer. 
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Royal  and  Bergles  (1978)  suggested  that  the  twisted-tape  Reynolds  number  (ReJ  replace  the  smooth  tube 
Reynolds  number  and  that  a fin-effect  correction  be  used.  The  fin-effect  correction  was  not  used  because 
later  work  by  Manglik  and  Bergles  (1987)  proved  it  to  be  negligible  for  snug  to  loose  fitting  tape  inserts. 
The  swirl  Reynolds  number  (Res)  is  based  on  the  swirl  velocity  (V^  and  the  hydraulic  diameter  of  the 
tube  with  the  twisted-tape  insert: 


(6) 


1 - 

7iDi 


For  most  of  the  data,  the  measured  condensation  heat-transfer  coefficient  decreases  for  decreasing 
qualities.  Apparently,  thin  liquid  films  and  high  vapor  velocities  at  the  entrance  of  the  condenser  provide 
for  high  heat-transfer  coefficients.  As  the  liquid  accumulates  on  the  tube  wall  for  decreasing  quality,  the 
heat-transfer  coefficient  diminishes.  The  decrease  in  the  heat-transfer  coefficient  from  a quality  of  0.4 
to  that  of  0.05  is  relatively  modest  compared  to  the  decrease  shown  by  the  predictive  correlations 
available  in  the  literature.  For  example,  all  of  the  predictions  by  correlations  from  the  literature  show 
at  least  a factor  of  two  increase  in  for  qualities  from  0.05  to  0.4.  Presumably,  the  modification  to  the 
smooth  tube  correlation  does  not  completely  account  for  the  swirl  enhancement  for  the  lower  qualities. 

The  Soliman  (1968),  the  Shah  (1979),  and  the  Traviss  et  al.  (1973)  correlations  were  nearly  always 
parallel,  intersecting  only  at  high  and  low  qualities.  Typically  these  correlations  were  within  15  % to  25  % 
of  one  another.  The  Soliman  (1968)  correlation  more  closely  predicted  the  measured  Nusselt  Numbers 
than  the  Shah  (1979),  and  the  Traviss  et  al.  (1973)  correlations  did.  The  Traviss  et  al.  (1973)  correlation 
consistently  gave  the  largest  overpredictions  for  xq  greater  than  0.4.  As  a group,  the  correlations 
typically  overpredict  the  measured  Nusselt  numbers  from  10%  to  100%  for  qualities  greater  than  0.4  with 
the  greatest  overpredictions  occurring  for  the  higher  qualities. 

Flow  Boiling 

Figures  13  through  21  show  the  measured,  local-averaged,  flow-boiling,  heat-transfer  coefficient  (hj  as 
a function  of  thermodynamic  quality.  The  data  are  compared  to  predictions  from  Three  different 
modified  smooth  tube  correlations.  The  smooth  tube  flow-boiling  correlations  of  Jung  et  al.  (1989), 
Kandlikar  (1990),  and  Gungor  and  Winterton  (1987)  were  modified  in  the  same  manner  as  the 
condensation  correlations  were  modified  to  account  for  the  effect  of  the  twisted  tape.  Namely,  the  Re,) 
was  replaced  with  Res  wherever  it  was  encountered  in  the  smooth  tube  correlations.  The  Jung  et  al. 
(1989)  correlation  was  the  only  one  valid  for  both  pure  components  and  mixtures.  The  pure  component 
correlations  of  Kandlikar  (1990),  and  Gungor  and  Winterton  (1987)  were  used  with  mixture  properties 
to  predict  the  mixture  heat  transfer.  A reasonable  fluid  factor  that  gave  a favorable  comparison  to  the 
data  was  used  in  the  Kandlidar  (1990)  correlation  for  the  mixture  predictions. 

The  shape  of  the  measured  versus  xq  plots  is  similar  for  the  various  fluids.  For  each  fluid,  the 
measured  h,  exhibits  two  distinct  regions  intersecting  at  approximately  xq  = 0.35.  The  bubbly  flow 
region  and  the  convective  flow  region  reside  for  approximately  Xq  < 0.35  and  for  Xq  > 0.35, 
respectively.  The  data  exhibits  a small  negative  slope  in  the  convective  region.  In  the  bubbly  flow 
region,  relatively  large  increases  in  the  heat-transfer  coefficient  occur  for  small  increases  in  the  quality. 
Speculation  on  the  basis  of  the  relationship  between  h,  and  Xq  in  the  convective  and  bubbly  flow  regions 
is  presented  in  the  following  four  paragraphs. 
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The  decrease  in  the  measured  h,  in  the  convective  region  can  be  attributed  to  increasing  partial  dryout  of 
the  tube  wall  with  increasing  quality.  The  twisted  tape  encourages  dryout  with  a 64%  greater  surface 
perimeter  than  that  of  an  empty  tube.  Consequently,  liquid  that  would  have  wetted  the  wall  of  the  empty 
tube  is  drawn/impelled  to  the  tape  by  surface-tension  and  vapor  acceleration  forces.  Dry  areas  on  the 
tube  wall  become  more  prevalent  as  less  liquid  is  available  to  wet  the  wall  and  larger  heat  fluxes  are 
available  to  evaporate  the  thin  liquid  films. 

The  promotion  of  partial  dryout  by  a twisted  tape  does  not  contradict  the  notion  that  twisted  tapes  can 
be  used  to  delay  burnout.  Bergles  et  al.  (1971)  conjecture  that  the  critical  heat  flux  condition  in  dispersed 
flow  is  delayed  as  a result  of  centrifugal  flow  forces  which  impel  small  liquid  droplets  to  the  outer  wall 
for  evaporation.  As  noted  by  Bergles  et  al.  (1971),  droplet  impinging  does  not  nullify  tape  wetting  even 
in  dispersed  flow.  Consequently,  the  twisted  tape  can  simultaneously  benefit  and  deter  heat  transfer  at 
high  qualities  by  delaying  burnout  and  inducing  partial  dryout,  respectively. 

The  correlations  from  the  literature  depict  a trend  that  is  contrary  to  that  of  the  measured  data.  The 
correlations  taken  from  the  literature  predict  an  increasing  heat-transfer  coefficient  with  quality  for  nearly 
the  entire  quality  range.  To  the  best  of  the  authors’  knowledge,  all  the  literature  correlations  were 
developed  with  data  primarily  from  electrical  resistance  heated  test  sections.  It  is  unlikely  that  any  partial 
dryout  data  was  included  in  these  correlating  data  sets.  Tube  dryout  for  an  electrically  heated  test  section 
would  present  a destructive  burnout  condition.  Consequently,  the  literature  correlations  present  a case 
where  the  walls  remain  wetted  as  the  quality  increases.  For  this  case,  increasing  vapor  shear  and 
decreasing  film  thickness  with  increasing  quality  act  to  enhance  the  heat  transfer. 

Figure  22,  a sample  plot  of  the  heat  flux  versus  wall  superheat  (ATj),  aids  in  explaining  the  versus 
quality  trend  in  the  bubbly  flow  region.  As  seen  in  Fig.  22,  the  heat  flux  is  relatively  small  but  rapidly 
increasing  in  the  low-quality  region.  Accordingly,  little  change  in  the  quality  occurs  in  the  bubbly  flow 
region  because  the  energy  input  is  small.  The  follows  the  rapidly  increasing  heat  flux  because  the  wall 
superheat  is  nearly  constant  in  this  region.  The  rapid  increase  in  the  heat  flux  is  a consequence  of  the 
counterflow  heat  exchange  between  the  refrigerant  and  the  secondary  fluid.  The  constant  superheat 
results  from  the  boiling  flow  pattern.  For  example,  little,  if  any,  bubble  nucleation  occurs  in  the  bubbly 
flow  region  because  of  the  small  heat  flux.  Instead,  tiny  bubbles  originate  from  the  expansion  valve. 
For  example,  Aaron  and  Domanski  (1989)  describe  the  flow  after  an  expansion  device  as  a misty  jet. 
The  ATS  remains  relatively  small  and  constant  in  the  bubbly  flow  region  because:  (1)  the  vapor  is 
thoroughly  interspersed  within  the  liquid  resulting  in  a near  equilibrium  state,  and  (2)  a negligible  amount 
of  superheat  is  required  to  maintain  and  grow  established  bubbles.  Liquid  superheat  occurs  for  qualities 
greater  than  0.35  because  the  vapor  has  coalesced  and  segregated  itself  from  the  liquid.  In  summary,  the 
rapidly  increasing  1^  in  the  low-quality  region  results  from  a rapidly  increasing  heat  flux  while  the 
refrigerant  liquid  at  the  wall  remains  near  saturation.  Consequently,  the  authors  anticipate  that  their 
boiling  correlation  would  be  valid  in  the  low-quality  range  for  refrigeration  cycles  where  bubble  mist  flow 
is  introduced  in  the  evaporator  in  the  absence  of  bubble  nucleation  at  the  wall. 

Considering  the  different  conditions  for  which  the  present  data  was  taken  and  the  literature  correlations 
were  derived,  it  is  not  surprising  that  the  predictions  and  the  measurements  differ.  The  predicted  pure 
fluid  Nusselt  numbers  obtained  from  the  modified  correlations  of  Jung  et  al.  (1989),  Kandlikar  (1990), 
and  Gungor  and  Winterton  (1987)  were  relatively  consistent  with  one  another.  For  the  most  part,  these 
correlations  over-predict  the  measured  Nusselt  numbers  for  convective  region  and  under-predict  those  for 
the  bubbly  flow  region.  It  was  not  uncommon  for  the  convective  region  to  be  over-predicted  by  100%. 
The  Nusselt  numbers  obtained  from  the  modified  correlations  typically  predicted  the  measured  Nusselt 
numbers  for  the  bubbly  flow  region  to  within  50%.  In  contrast  to  predictions  for  the  pure  fluids,  the 
azeotropic  predictions  from  the  Kandlikar  (1990)  correlation  were  as  much  as  30%  lower  than  those 
obtained  from  the  Jung  et  al.  (1989),  and  the  Gungor  and  Winterton  (1987)  correlations  which  remained 
consistent  with  each  other. 
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CORRELATIONS 

From  the  previous  section  it  is  obvious  that  the  modified  existing  correlations  do  not  reliably  predict  the 
present  data.  The  following  section  presents  three  correlations  that  were  derived  from  the  present  Nusselt 
number  measurements.  The  three  correlations  quantify  the  local  Nusselt  numbers  for  convective  boiling 
and  condensation  of  pure/azeotropic  and  zeotropic  fluids  within  a tube  with  a twisted-tape  insert.  The 
fluid  properties  and  heat  transfer  conditions  used  in  the  correlations  were  evaluated  locally. 

Manglik  and  Bergles  (1992)  suggest  that  swirl  flow  effects  may  be  correlated  by  the  swirl  parameter 
(Sw): 


Sw 


Res 

v/y 


(7) 


The  1 /v/y  parameter  accounts  for  convective  inertia  effects  that  are  ignored  by  1/y  which  has  been 
traditionally  used  to  correlate  swirl  flow  data.  Manglik  and  Bergles  (1992)  show  that  laminar  flow  heat 
transfer  for  three  different  twist  ratios  are  correlated  well  with  the  swirl  parameter.  Agrawal  et  al.  (1986) 
also  correlated  R12  swirl  flow-boiling  heat-transfer  coefficients  to  y05219  for  ReD  approximately  in  the 

range  of  7000  to  14000.  Consequently,  it  is  inferred  that  the  lly/y  parameter  is  suitable  for  correlating 
flow  boiling  and  convective  condensation  data  for  low  and  somewhat  high  Reynolds  numbers. 

The  Law  of  Corresponding  States  philosophy  presented  by  Cooper  (1984)  was  followed  by  including  the 
reduced  pressure  (Pr),  -log10Pr,  and  the  molecular  weight  (M)  in  the  correlation.  Cooper  (1984)  suggests 
that  the  fluid  properties  that  govern  nucleate  pool  boiling  can  be  well  represented  by  a product  of  these 
variables  to  various  powers.  The  above  terms  and  several  others  were  used  to  correlate  our  data  for  all 
conditions  of  two-phase  flow. 

The  Nusselt  number  correlations  were  derived  from  phase  change  data  taken  in  a 9.6  mm  diameter  tube 
with  a twisted-tape  of  twist  ratio  4.15,  and  for  the  parameter  ranges  given  in  Table  3.  The  heat  transfer 
measurements  were  fitted  to  three  separate  correlations.  The  convective  condensation  Nusselt  numbers 
for  single  component,  azeotropic,  and  zeotropic  fluids  were  fitted  to  a single  equation.  Two  equations 
were  necessary  to  correlate  the  flow  boiling  Nusselt  numbers  of  refrigerants  and  refrigerant  mixtures. 
The  data  were  fitted  to  fundamental  dimensionless  parameters  in  hopes  that  the  correlations  could  be 
extrapolated  to  other  fluids,  tube  diameters,  and  twist  ratios.  There  is  no  way  to  know  a priori  how  well 
the  equations  will  perform  when  they  are  extrapolated  beyond  the  conditions  of  the  data  base. 

There  are  several  details  that  are  common  to  all  three  of  local  Nusselt  number  correlations  that  follow. 
First,  all  of  the  parameters  were  evaluated  locally.  The  range  of  each  parameter  as  correlated  is  given 
in  Table  3.  Each  parameter  is  raised  to  a nonconstant  exponent  which  is  a quadratic  function  of  quality. 
The  quality  dependent  exponents  enable  the  single  correlation  to  predict  most  of  the  two-phase  region. 
For  example,  as  will  be  seen  in  eqn  8,  the  exponent  on  the  boiling  number  (Bo)  is  relatively  large  in  the 
low  quality  region  and  decreases  for  the  high  quality  region.  Finally,  all  Nusselt  number  measurements 
with  a relative  expanded  uncertainty  greater  than  12.5%  where  omitted  from  the  regression.  The 
omissions  resulted  in  a loss  of  approximately  1 % of  the  data. 

Table  3 shows  that  the  boiling  data  was  taken  for  transition  and  turbulent  all-liquid  Reynolds  numbers 
(ReD).  The  all-liquid  Reynolds  were  evaluating  using  the  local  liquid  properties  of  the  fluid.  The  table 
also  shows  that  the  condensation  data  were  taken  for  laminar  and  turbulent  all-liquid  Reynolds  numbers. 


Flow-boiling 

The  1401  locally  measured  flow-boiling  Nusselt  numbers  for  refrigerants  R12,  R22,  R152a,  R134a, 
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R290,  R290/R134a,  and  R134a/R600a  were  correlated  to  the  following  equation: 


NU_  = 1.356  * Sw 


0.993-1. 181x„+Q . 899xr 


pr  1. 108-2. 366xq+1.451xq2 
p -2 . 383+5 .255xq-l .791xq2 
-3 . 195+6 . 66  8x~ 


(8) 


• (-log10Pr) 


^ 1 . 073-2 .679xct+1.443x2 

BO  q q 


Figure  23  shows  that  eqn  8 correlates  95  % of  the  pure  component  and  azeotropic  flow-boiling  Nusselt 
numbers  to  within  approximately  ± 25  % . The  mean  of  the  measurements  has  an  expanded  uncertainty 
of  ± 5 % . Figures  24  through  30  show  the  measured  versus  predicted  Nusselt  numbers  for  each  fluid 
along  with  the  95%  confidence  interval  for  the  mean  of  the  measurements.  Only  random  trends  were 
observed  in  the  residual  plots  against  each  of  the  parameters  of  eqn.  8.  Equation  8 is  plotted  as  a solid 
line  on  figures  13  through  19. 

Figure  31  compares  eqn.  8 to  measured  flow-boiling  heat-transfer  coefficients  from  the  literature. 
Agrawal  et  al.  (1986)  present  local  R12  flow-boiling  heat-transfer  coefficients  for  a tube  with  a twisted 
tape  with  a twist  ratio  of  5.58.  Equation  8 predicts  the  Agrawal  et  al.  (1986)  data  to  within  15%  for 
qualities  between  30%  and  40%.  Larger  differences  between  Equation  8 and  the  Agrawal  et  al.  (1986) 
data  were  evident  for  higher  heat  flux  conditions.  Figure  31  also  compares  the  locally  measured  water 
flow-boiling  heat-transfer  coefficients  by  Blatt  and  Adt  (1963)  to  eqn.  8 for  a twist  ratio  of  5.  A precise 
assessment  of  eqn.  8 using  the  Blatt  and  Adt  (1963)  data  is  difficult  due  to  the  large  scatter  in  the  data. 
Also,  Blatt  and  Adt  (1963)  did  not  give  the  test  pressure  in  their  manuscript.  As  a result,  eqn.  8 is 
plotted  using  two  likely  test  pressures.  The  most  that  can  be  said  is  that  the  correlation  passes  through 
the  vicinity  of  the  data. 


Schluender  (1982)  developed  a nucleate  pool  boiling  model  for  binary  zeotropic  mixtures.  His  model 
includes  the  effects  of  Bo,  the  difference  in  the  vapor  and  liquid  mole  fraction  of  the  more  volatile 
component  (x^  - x,),  and  the  difference  in  the  saturation  temperature  of  the  two  pure  components  (TLV  - 
Tmv.)  Thome  (1989)  criticizes  the  use  of  TLV  - Tmv  and  states  that  the  mixture’s  boiling  range  would  be 
a more  appropriate  parameter  for  the  model.  Accordingly,  the  present  correlations  contain  the  mole 
fraction  difference,  the  Boiling  number,  and  a dimensionless  mixture  temperature  difference  (0): 


0 = 


(9) 


where  Tb  and  Td  are  the  local  bubble  point  and  dew  point  temperatures  of  the  mixture,  respectively. 
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The  935  zeotropic  flow-boiling  Nusselt  numbers  for  R32/R152a  and  R32/R134a  were  fitted  to: 


Num 

Nup 


= 2.525  * (xv  - x^-°-58+5-67xq-2-825x. 

. 0-2.793xq  # p 1. 204-3. 335xq+1.946xq2 
* BO  0,338+’353xq  . .839 


(10) 


where  Niip  was  evaluated  using  eqn.  8 and  the  local  properties  of  the  zeotropic  mixture. 

Figure  32  shows  that  eqn  10  correlates  95%  of  the  zeotropic  flow-boiling  Nusselt  numbers  to  within 
approximately  ± 20% . The  mean  of  the  measurements  has  an  expanded  uncertainty  of  ± 4% . Figures 
33  and  34  show  the  measured  versus  predicted  Nusselt  numbers  for  each  fluid  along  with  the  95% 
confidence  interval  for  the  mean  of  the  measurements.  Only  random  trends  were  observed  in  the  residual 
plots  against  each  of  the  parameters  of  eqn  10.  Equation  10  is  plotted  as  a solid  line  on  figures  20  and 
21. 


Convective  condensation 

The  2253  locally  measured  convective  condensation  Nusselt  numbers  for  refrigerants  R12,  R22,  R152a, 
R134a,  R290,  R290/R134a,  R134a/R600a,  R32/R134a,  and  R32/R152a  were  correlated  to  the  following 
equation: 


Nu  = 0 . 00136  Sw 


0 .613+0 .647x, 


pr0.877.  p -1. 735+2. 362xq 


• (-log10Pr) 

• Ja~°-528 


-2 . 815+4 . 197x, 


(ID 


Figure  35  shows  that  eqn  11  correlates  95%  of  all  the  condensation  Nusselt  numbers  to  within 
approximately  ± 20%.  The  single  equation  was  found  to  acceptably  predict  single  component  and 
zeotropic  Nusselt  numbers  equally  well.  This  suggests  that  mass  transfer  has  a neglible  effect  on 
condenstion  heat  transfer  with  twisted  tape-inserts.  The  mean  of  the  measurements  has  an  expanded 
uncertainty  of  ± 3 % . 

Figures  36  through  45  show  the  measured  versus  predicted  Nusselt  numbers  for  each  fluid  along  with  the 
95%  confidence  interval  for  the  mean  of  the  measurements.  Only  random  trends  were  observed  in  the 
residual  plots  against  each  of  the  parameters  of  eqn  1 1 . Equation  1 1 is  plotted  as  a solid  line  on  figures 
4 through  12. 

CONCLUSIONS 

The  Nusselt  numbers  for  flow  boiling  and  convective  condensation  with  a twisted-tape  insert  were 
measured  and  correlated  for  several  single  component,  azeotropic,  and  zeotropic  refrigerants.  Three 
different  correlations  were  developed  to  predict  the  flow  boiling  and  convective  condensation  conditions 
for  azeotropic  and  zeotropic  refrigerants.  The  Nusselt  numbers  were  correlated  with  dimensionless 
variables  raised  to  nonconstant  exponents.  The  nonconstant  exponents  enabled  each  single  correlation 
to  capture  the  variation  of  the  relative  importance  of  each  parameter  with  quality.  The  flow-boiling 
correlations  are  valid  for  transition  and  turbulent  all-liquid  Reynolds  numbers.  The  convective 
condensation  correlation  is  valid  for  laminar  and  turbulent  all-liquid  Reynolds  numbers,  and  for  pure. 
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azeotropic  and  zeotropic  refrigerants. 

The  condensation  measurements  were  compared  to  several  smooth  tube  correlations  that  were  modified 
to  account  for  the  effect  of  swirl  flow.  For  the  low  quality  region,  the  rate  of  decrease  in  the  twisted-tape 
condensation-Nusselt  number  was  significantly  less  than  that  predicted  with  the  modified  correlations. 
Also,  the  magnitude  of  the  measured  condensation-Nusselt  numbers  were  significantly  greater  than  the 
predictions  for  the  low  quality  region.  Presumably,  the  modification  to  the  smooth  tube  correlation  does 
not  completely  account  for  the  swirl  enhancement  for  the  lower  qualities. 

The  flow-boiling  measurements  were  also  compared  to  several  smooth  tube  correlations  that  were 
modified  to  account  for  the  effect  of  swirl  flow.  For  the  most  part,  the  modified  correlations  significantly 
over  predicted  the  measured  boiling-Nusselt  numbers  in  the  convective  region.  It  was  conjectured  that 
the  twisted  tape  encouraged  partial  dryout  of  the  tube  wall.  The  partial  dryout  caused  the  Nusselt  number 
to  decrease  with  increasing  quality.  The  modified  correlations  predicted  an  opposite  trend  with  quality. 
Possibly,  the  fact  that  the  modified  correlations  were  derived  from  data  obtained  from  electrical  resistance 
heated  rigs  prohibits  the  prediction  of  the  dryout  effect.  Also,  the  very  low  quality  heat  transfer 
measurements  were  not  predicted  very  well  by  the  modified  correlations.  An  expansion  valve  supplied 
a bubbly-mist  flow  to  the  entrance  of  the  evaporator.  Consequently,  no  nucleation  of  bubbles  occurred 
at  the  wall.  This  is  in  contrast  to  electrically  heated  test  rigs  where  high  superheats  and  bubble  nucleation 
at  low  quality  are  common. 
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Table  1 R32  Mass  Compositions  for  Zeotropic  Mixtures 


% Mass 

Composition  of 
R32 

R32/134a 

R32/152a 

15 

X 

17 

X 

19 

X 

20 

X 

21 

X 

X 

22 

X 

23 

X 

24 

X 

X 

25 

X 

26 

X 

X 

29 

X 

30 

X 

X 

31 

X 

X 

32 

X 

33 

X 

34 

X 

X 

35 

X 

X 

36 

X 

X 

37 

X 

X 

38 

X 

39 

X 

X 

40 

X 

X 

41 

X 

X 

1 42 

X 

43 

X 
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Table  2 Average  Expanded  Relative  Measurement  Uncertainties  (%) 


qe 

qc 

Re 

Pr 

Bo 

UA 

X 

m 

Tur 

bine 

m 

Cori 

olis 

h* 

q"e 

q"c 

Ja  I 

4 

4 

3 

0.5 

8.0 

7 

7 

10 

3 

17 

5 

8 

3 
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Table  3 Parameter  Range  for  Correlations 


parameter 

flow  boiling 
pure  and 
azeotropic 

flow 

condensation 
pure  and 
azeotropic 

flow 

boiling 

zeotropic 

mixture 

flow 

condensation 

zeotropic 

mixture 

Re 

1800  - 15000 

4000-32600 

2600- 

12000 

5200  - 18500 

Sw 

1000-8430 

2248-18300 

1460- 

6750 

3000  - 10400 

Pr 

0.035-0.2 

0.14-0.45 

0.032  - 
0.14 

0.15  - 0.34 

Pr 

2.4 -4.3 

2.4 -3.6 

2.2  - 
3.8 

2.2  - 3.3 

Bo 

2.8  x ia5 

5.i  x ia4 

n/a 

2.5  x ia5 

4.8  x 104 

n/a 

Ja 

n/a 

0.008-0.1 

n/a 

0.007-0.05 

M 

44-121 

44-121 

61-92 

61-92 

Xq 

0.11  -0.98 

0.023  - 0.9 

0.12  - 1 

0.016-0.9 

Xy-Xi 

n/a 

n/a 

0.12  - 
0.26 

0.1  -0.2 

ATb 

n/a 

n/a 

0.14  - 
0.25 

0.1  - 0.2 
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HC 

AAA 

■ > (T)"  . 

_zfx_  Three-way 
''■f'  valve 

_v7_j<p\  Thermostatic 
~ valve,  electric 

V Thermostatic 
^ ^ valve 

I^|  Hand  valve 

0 Filter,  dryer 
(refrigerant) 

© Sight  glass 
Filter  (water) 


[WV 

WV 

^ \ 

HC 

HV 

C - Condenser 
Co  - Compressor 
E - Evaporator 
FM  - Volume  flow  meter 
MM  - Mass  flow  meter 


He  - Heater,  constant 
HV  - Heater,  variable 
I - Inverter 

IX  - Internal  heat  exchanger 
M - Motor 


P - Pump 

S - Torque  & speed  meter 
W/G  - Water/glycol  tank 
X - External  heat  exchanger 


Fig.  1 Schematic  of  the  breadboard  refrigeration  cycle  (Pannock 
and  Didion,  1991) 
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Fig.  2 Schematic  of  evaporator  or  condenser  with  thermocouple 
locations 
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Fig.  3 Refrigerant  and  secondary  fluid  temperature  profile  of 
R22  evaporator  and  condenser 
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Fig.  4 Convective  condensation,  twisted- tape  Nusselt  number  as 
a function  of  thermodynamic  quality  for  R12 
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Fig.  5 Convective  condensation,  twisted- tape  Nusselt  number  as 
a function  of  thermodynamic  quality  for  R134a 
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Fig.  6 Convective  condensation,  twisted- tape  Nusselt  number  as 
a function  of  thermodynamic  quality  for  R22 
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Fig.  7 Convective  condensation,  twisted-tape  Nusselt  number  as 
a function  of  thermodynamic  quality  for  R290 
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Fig.  8 Convective  condensation,  twisted- tape  Nusselt  number  as 
a function  of  thermodynamic  quality  for  R134a/R600a 
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Fig.  9 Convective  condensation,  twisted- tape  Nusselt  number  as 
a function  of  thermodynamic  quality  for  R290/R134a 


28 


600 

500 

400 

Nu  300 

200 

100 

0 

0 0.1  0.2  0.3  0.4  0.5  0.6  0.7  0.8  0.9  1.0 


Fig.  10  Convective  condensation,  twisted-tape  Nusselt  number  as 
a function  of  thermodynamic  quality  for  R152a 
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Fig.  11  Convective  condensation,  twisted- tape  Nusselt  number  as 
a function  of  thermodynamic  quality  for  R32/R134a 
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Fig.  12  Convective  condensation,  twisted- tape  Nusselt  number  as 
a function  of  thermodynamic  quality  for  R32/R152a 
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Fig.  13  Flow  boiling,  twisted- tape  Nusselt  number  as  a function 
of  thermodynamic  quality  for  R12 
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Fig.  14  Flow  boiling,  twisted- tape  Nusselt  number  as  a function 
of  thermodynamic  quality  for  R134a 
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Fig.  15 


Flow  boiling,  twisted-tape  Nusselt  number  as  a function 
o±  thermodynamic  quality  for  R22 
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Fig.  16  Flow  boiling,  twisted- tape  Nusselt  number  as  a function 
of  thermodynamic  quality  for  R290 
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Fig.  17  Flow  boiling,  twisted-tape  Nusselt  number  as  a function 
of  thermodynamic  quality  for  R134a/R600a 
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Fig.  18  Flow  boiling,  twisted- tape  Nusselt  number  as  a function 
of  thermodynamic  quality  for  R290/R134a 
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Fig.  19  Flow  boiling,  twisted-tape  Nusselt  number  as  a function 
of  thermodynamic  quality  for  R152a 
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Fig.  20  Flow  boiling,  twisted-tape  Nusselt  number  as  a function 
of  thermodynamic  quality  for  R32/R134a 
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Fig.  21  Flow  boiling,  twisted-tape  Nusselt  number  as  a function 
of  thermodynamic  quality  for  R32/R152a 
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Fig.  22  Sample  plot  of  the  heat  flux  versus  wall  superheat 
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Fig.  23  Correlation  of  twisted- tape  flow  boiling  of  single 
components  and  azeotropic  mixture  Nusselt  numbers 
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Fig.  24  Correlation  of  twisted-tape , R12,  flow-boiling  Nusselt 
numbers 
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Fig.  25  Correlation  of  twisted-tape , R134a,  flow-boiling  Nusselt 
numbers 
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Fig.  26  Correlation  of  twisted-tape , R22,  flow-boiling  Nusselt 
numbers 
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Fig.  27  Correlation  of  twisted-tape , R290,  flow-boiling  Nusselt 
numbers 
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Fig.  28  Correlation  of  twisted-tape , R134a/R600a  (0.81/0.19), 

flow-boiling  Nusselt  numbers 
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Fig.  29  Correlation  of  twisted-tape,  R290/R134a  (0.44/0.56), 

flow-boiling  Nusselt  numbers 
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Fig.  30  Correlation  of  twisted-tape , R152a,  flow-boiling  Nusselt 
numbers 
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Fig.  31  Comparison  of  twisted-tape  flow  boiling  correlation  of 
data  available  in  the  literature 
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Fig.  32  Correlation  of  twisted- tape  flow  boiling  of 
mixture  Nusselt  numbers 
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Fig.  33  Correlation  of  R32/R134a,  twisted-tape , flow 
Nusselt  numbers 
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Fig.  34  Correlation  of  R32/R152a,  twisted-tape , flow-boiling 
Nusselt  numbers 
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Fig.  35  Correlation  of  twisted-tape,  convective  condensation  of 
sin^ls  component  and  azeotropic  mixture  Nusselt  numbers 
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Fig.  36  Correlation  of  twisted- tape , convective 
zeotropic  mixture  Nusselt  numbers 
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Fig.  37  Correlation  of  R32/R134a,  twisted-tape , 
condensation  Nusselt  numbers 
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Fig.  38  Correlation  of  R32/R152a,  twisted-tape , 
condensation  Nusselt  numbers 
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Fig.  39  Correlation  of  twisted-tape,  R12,  convective  condensation 
Nusselt  numbers 
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Fig.  40  Correlation  of  twisted-tape,  R134a,  convective 
condensation  Nusselt  numbers 
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Fig.  41  Correlation  of  twisted- tape,  R22,  convective  condensation 

Nusselt  numbers 
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Fig.  42  Correlation  of  twisted-tape , R152a,  convective 

condensation  Nusselt  numbers 
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Fig.  43  Correlation  of  twisted- tape , R290,  convective 

condensation  Nusselt  numbers 
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Fig.  44  Correlation  of  twisted-tape , R134a/R600a,  convective 

condensation  Nusselt  numbers 
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Fig.  45  Correlation  of  twisted-tape,  R290/R134a,  convective 
condensation  Nusselt  numbers 
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APPENDIX  A1 


Flow  boiling  of  pure  components  and 
azeotropes  with  a twisted  tape 

(file:  combe.dat) 
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APPENDIX  A2 


Convective  condensation  of  pure 
and  azeotropes  with  a twisted  tape 
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Flow  boiling  of  zeotropes 
with  a twisted  tape 
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